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SHIP-TRACK  MODELS  BASED  ON  POISSON  DISTRIBUTED 
PORT-DEPARTURE  TIMES 


1.  INTRODUCTION 

A  number  of  applications  require  a  statistical  description  of  the  locations  of  ships  in  a  region,  e.g., 
ambient  noise  models  and  probability-based  ship  tracking  algorithms.  For  many  of  these  applications,  it 
suffices  to  specify  the  mean  number  of  ships  on  a  grid  of  elemental  longitude-latitude  areas.  Such 
descriptions,  often  referred  to  as  shipping  distributions,  have  been  obtained  from  both  direct 
measurements  of  individual  ship  positions  and  from  shipping  models  based  on  port-departure  times  and 
routing  data.  A  well-known  example  of  the  latter  is  the  Historical  Temporal  Shipping  (HITS)  model  [1,2]. 
This  model  specifies  shipping  distributions  for  different  classes  of  commercial  ships  and  different  time 
periods  over  a  large  portion  of  the  world's  oceans.  The  resolution  of  these  shipping  distributions,  i.e.,  the 
dimensions  of  the  longitude-latitude  grid,  is  1°  x  1°. 

Many  applications,  however,  require  a  description  of  the  tracks  of  the  individual  ships  over  some  time 
interval  (t0,t0  +  7']  rather  than  simply  the  mean  number  of  ships  in  elemental  areas.  For  some  of  these 
applications,  the  duration  of  the  required  observation  period  T  may  exceed  the  time  interval  over  which 
many  of  the  ships  are  within  the  region  of  interest.  For  these  applications,  the  ship-track  model  must  be 
capable  of  introducing  new  ships  into  the  region  in  a  manner  consistent  with  the  underlying  shipping 
distribution. 

The  requirements  of  the  ship-track  model  also  depend  on  the  types  of  ships  that  are  relevant  to  the 
application.  For  most  commercial  shipping,  the  ships  simply  transit  from  one  port  to  another  along  routes 
that  constitute  a  more  or  less  well-defined  shipping  lane.  For  this  type  of  shipping,  it  may  suffice  to 
provide  simple  smooth  approximations  to  the  actual  routes  traveled  from  the  departure  to  the  destination 
port.  On  the  other  hand,  shipping  traffic  such  as  fishing,  recreational  and  military  vessels  can  travel  along 
complicated  routes  as  dictated  by  the  specific  mission  of  the  vessel.  For  this  type  of  shipping,  more 
complicated  descriptions  of  the  routes  may  be  required.  In  either  case,  the  accuracy  requirements  of  the 
ship  tracks  will  depend  on  the  specific  application  of  the  model. 

This  report  presents  two  models,  a  deterministic  and  a  stochastic  model,  each  of  which  describe  the 
tracks  of  ships  en  route  in  a  region  during  an  arbitrary  time  interval  (tQ,t0  +T] .  Both  models  consist  of  a 
track  function  that  describes  the  tracks  of  the  individual  ships  and  a  probability  law  on  the  total  number  of 
ships  en  route  during  (tQ,t0  +  7  ] ,  the  positions  of  those  ships  at  the  initial  time  t0 ,  and  their  nominal 
speeds.  The  probability  law  is  obtained  under  the  assumption  that  the  times  at  which  the  ships  depart  each 
port  are  Poisson  distributed  with  a  time-varying  departure  rate  and  that  the  ship  speeds  and  the  routes  that 
the  ships  travel  are  statistically  independent.  Under  this  assumption,  it  is  shown  that  the  ship-track 
parameters  are  distributed  as  a  Poisson  process  with  a  time-dependent  rate  function.  The  rate  function  that 
specifies  this  process  differs  for  the  two  ship-track  models. 

The  two  ship-track  models  provide  alternate  descriptions  of  the  ship’s  tracks.  In  the  deterministic 
model,  the  ship  tracks  are  deterministic  functions  derived  from  a  probability  density  on  ship  positions  in 
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manner  analogous  to  that  of  the  HITS  model  (see  Ref.  2).  As  such,  the  individual  ship  tracks  can  neither 
double  back  on  themselves  nor  intersect  any  other  ship  tracks.  This  may  be  an  acceptable  constraint  for 
ship  traffic  that  simply  transits  from  one  port  to  another,  provided  that  detailed  representations  of  the 
individual  tracks  are  not  required.  In  the  stochastic  model,  the  ship  tracks  are  obtained  as  realizations  of  a 
Markov  process  without  any  constraints  on  the  route  crossings.  Accordingly,  this  model  is  suitable  for 
shipping  where  the  tracks  are  determined  by  operations  more  complicated  than  port-to-port  transiting. 
Furthermore,  the  stochastic  nature  of  the  tracks  may  be  more  realistic  than  the  smooth  track 
approximations  of  the  deterministic  model. 

This  report  is  organized  as  follows:  Sections  2  and  3  provide  the  background  for  the  model 
definitions.  In  Section  2,  we  present  the  definitions  and  assumptions  that  are  used  to  describe  the 
kinematics  of  the  shipping  for  each  ordered  pair  of  ports  that  support  shipping  in  the  region.  Specifically, 
for  each  port  pair,  we  define  the  ship-track  function  in  terms  of  a  route  function  and  a  motion  function.  In 
general,  these  functions  depend  on  a  track  parameter  that  specifies  the  distance  that  the  ship  has  traveled 
at  the  initial  time  t0 ,  the  specific  route  that  the  ship  travels  between  the  departure  and  the  destination  port 
and  the  ships  nominal  speed  along  that  route.  To  account  for  those  ships  that  are  not  en  route  at  time  t0 , 
but  depart  the  port  during  the  interval  (t0 ,  t0  +  7'j ,  we  allow  these  parameters  to  take  on  “virtual”  values. 
In  Section  3,  we  state  the  assumptions  on  the  port-departure  probability  law  and  describe  the  probability 
law  on  the  track  parameter  that  results  from  these  assumptions.  It  is  seen  that  the  track-parameter 
probability  law  can  be  viewed  as  the  composition  of  two  Poisson  process,  one  representing  the  ships  that 
are  present  at  time  t0  and  one  representing  the  ships  that  enter  the  region  in  the  interval  (t0 ,t0  +  T\ . 

The  ship-track  models  themselves  are  described  in  Sections  4  and  5.  For  each  of  these  models,  we 
first  specify  the  route  set  probability  law  that  forms  the  basis  of  the  model.  We  then  specify  the  track 
function  and  the  track-parameter  probability  law  that  results  from  the  route  set  probability  law.  For  both 
models,  the  track  parameter  can  be  expressed  in  terms  of  the  ship  position  at  the  initial  time  t0  and  the 
nominal  speed.  This  leads  to  the  notion  of  the  shipping  density  from  which  the  shipping  distribution  can 
be  derived  by  integrating  over  elemental  areas.  The  derivations  of  these  results  are  presented  in 
Appendixes  A  and  B.  For  reference  purposes  we  have  also  provided  a  definition  of  a  multi-dimensional 
Poisson  process  in  Appendix  A.  In  Appendix  C,  we  present  a  simplified  version  of  the  deterministic 
model.  Finally,  the  results  of  the  report  are  summarized  in  Section  6.  An  example  of  an  ambient  noise 
model  that  draws  on  a  deterministic  ship-track  model  is  described  in  Refs.  3  and  4.  Reference  5  presents 
an  ambient  noise  application  to  the  region  around  San  Diego,  California. 

2.  KINEMATICS 

A  realization  of  ship  tracks  must  represent  not  only  those  ships  that  are  present  in  the  region  during 
the  entire  time  interval  of  interest,  but  also  those  ships  that  either  enter  the  region  or  exit  the  region  during 
that  interval.  The  latter  do  so  by  either  departing  or  arriving  at  a  port  or  by  crossing  the  boundary  of  the 
region.  Figure  1  illustrates  this  process.  The  region  of  interest  is  bounded  by  the  coastline  segments 
indicated  by  the  heavy  black  lines  and  by  the  dotted  line  segments  across  the  access  areas  to  the  region. 
The  ports,  labeled  P\  through  P9 ,  are  classified  as  either  “real  ports”  or  “pseudo  ports.”  The  real  ports  are 
those  that  are  physically  located  within  the  region.  The  pseudo-ports,  located  along  the  access  boundaries 
of  the  region,  represent  the  shipping  from  the  aggregate  of  the  ports  located  outside  the  region  that 
provide  traffic  to  the  region.  The  ports  R6  and  P9  are  pseudo  ports;  the  remaining  ports  are  real  ports. 
The  ship  tracks  are  indicated  by  the  line  segments.  The  dots  on  the  segments  show  the  ship  positions  at 
time  tQ ;  the  arrowheads  show  the  positions  at  time  t0+T  .  The  ship  tracks  fall  into  two  groups,  ships  that 
are  en  route  at  time  t0  and  ships  that  are  in  port  at  time  t0  but  depart  the  port  during  the  interval 
(t0,t0  +T ] .  Ship  tracks  in  the  first  group  have  a  dot  at  the  beginning  of  the  track;  ship  tracks  in  the 
second  group  do  not  have  a  dot  at  the  beginning  of  the  track.  For  tracks  in  the  first  group,  there  are  a 
number  of  possibilities:  (a)  the  ship  remains  en  route  for  the  entire  time  interval  (indicated  by  a  track  with 
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both  a  dot  and  an  arrowhead);  (b)  the  ship  arrives  at  a  port  during  the  time  interval  and  remains  there  for 
the  rest  of  the  time  interval  (indicated  by  a  track  with  a  dot  but  no  arrowhead);  (c)  the  ship  arrives  at  a 
port  and  then  departs  from  that  port  during  the  time  interval  (indicated  by  the  dashed  track  at  P2  );  and  (d) 
the  ship  arrives  at  and  then  departs  from  more  than  one  port  during  the  interval.  The  same  possibilities 
exist  for  the  ships  that  are  at  a  port  at  the  beginning  of  the  time  interval  once  those  ships  have  departed  the 
port. 


Fig.  1  —  A  possible  set  of  ship  tracks  for  a  region 

The  ship-track  model  presented  here  does  not  describe  the  tracks  of  ships  that  pass  through  a  port 
during  the  time  interval  of  interest.  Instead,  it  represents  such  tracks  as  two  tracks  —  one  associated  with 
a  ship  that  arrives  at  the  port  and  remains  there  for  the  rest  of  the  time  interval  and  one  associated  with  a 
ship  that  is  at  the  port  at  the  beginning  of  the  interval  and  then  departs  at  some  time  during  the  interval. 
Thus,  the  track  of  the  ship  passing  through  P2  in  Fig.  1  is  represented  as  the  track  of  one  ship  that  enters 
P2  plus  the  track  of  another  ship  that  departs  from  P2  .  This  approximation  simplifies  the  model 
considerably  since  it  obviates  the  need  to  relate  a  ship’s  departure  time  to  its  arrival  time.  Furthermore,  it 
allows  the  shipping  between  any  pair  of  ports  to  be  considered  separately  from  the  shipping  between  any 
other  pair  of  ports.  This  approximation  is  reasonable  if  the  duration  of  the  time  interval  is  short  compared 
to  the  mean  time  that  ships  spend  in  a  port  since,  in  this  case,  the  number  of  tracks  that  pass  through  the 
port  in  the  time  interval  will  be  negligible.  It  is  also  be  reasonable  for  longer  time  intervals  if  it  is  not 
necessary  to  preserve  the  identity  of  the  individual  ships  that  pass  through  the  ports. 

As  a  result  of  the  independent  port-pair  approximation,  the  regional  shipping  model  can  be  obtained 
as  the  composition  of  ship-track  models  for  each  ordered  pair  of  ports  ( Pn ,  Pm )  that  support  shipping 
within  the  region.  The  first  port  in  the  port-pair  is  the  departure  port ;  the  last  port  is  the  destination  port. 
We  also  include  the  port-pair  (Pn,Pn)  since  for  certain  kinds  of  shipping  (e.g.,  fishing  and  recreational) 
the  destination  port  and  the  departure  port  can  be  the  same.  Furthermore,  we  include  both  the  port-pairs 
(Pn ,  Pm )  and  (Pm ,  Pn ) ,  since  the  traffic  from  Pn  to  Pm  is  not  necessarily  the  same  as  the  traffic  from  Pm 
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to  Pn  .  Thus,  if  there  are  Np  ports  for  the  region,  including  pseudo  ports,  there  are  as  many  as  N~  port- 
pair  models  for  the  regional  model. 

To  describe  the  kinematics  for  an  individual  port  pair,  it  is  first  necessary  to  specify  the  routes  that  the 
ships  travel  and  the  motion  of  those  ships  along  those  routes.  These  quantities  determine  a  ship-track 
function  that  specifies  the  tracks  of  the  individual  ships  during  the  period  of  interest. 

2.1  The  Route  Function 

The  routes  that  the  ships  travel  are  specified  by  a  route  function  that  describes  the  set  of  all  possible 
routes  from  Pn  to  Pm  by  associating  each  possible  route  with  a  route  parameter  co  e  Q  and  expressing 
that  route  as  a  function  of  the  distance  x  measured  along  that  route  from  the  departure  port.  The  route 
function  has  the  form 


where  R^(x;oj)  and  /^(xco)  are  the  coordinate  functions  for  the  longitude  coordinate  A  and  the  latitude 
coordinate  4* ,  respectively.  For  each  value  of  the  route  parameter  co ,  we  denote  the  length  of  the  route  by 
L( co)  and  assume  that  as  x  increases  from  0  to  L( co) ,  the  route  function  traces  out  the  route  co  from  Pn 
to  Pm  as  a  continuous  path  in  the  set  A  e  (— tc,tt]  and  <|)  e  (— tx /  2,  tx /  2) .  This  set  represents  all  positions  on 
the  Earth,  except  for  the  north  and  the  south  poles,  with  the  convention  that  negative  values  of  A 
correspond  to  longitudes  east  of  the  Greenwich  meridian  and  negative  values  of  4*  correspond  to  latitudes 
south  of  the  equator. 

To  account  for  ships  that  depart  from  Pn  after  time  t0 ,  we  allow  the  distance  traveled  x  to  be 
negative  and  define  the  route  function  by 

Ri  (x;co)  =  A-rc 

(2) 

(x;  co)  =  <()  =  co 

With  this  definition,  route  segments  determined  for  negative  values  of  x  do  not  describe  actual  routes 
since  the  coordinates  (A., 4>)  take  values  in  the  set  SJ?V  =  (-co,-7t)xfi .  We  refer  to  this  set  as  the  virtual 
ship  coordinates  and  use  it  to  implicitly  determine  the  departure  times  of  the  ships  that  leave  Pn  during 
the  time  interval  (t0,t0  +  T],  We  refer  to  the  set  iRr  =  (-7i,7t]x(-7t/2,7t/2)  as  the  real  ship  coordinates 
and  to  the  union  of  both  sets  =  '.H v  usJfr  as  the  extended  ship  coordinates. 

In  the  sequel,  it  is  convenient  to  define  the  route  function  and  the  track  function  in  an  auxiliary 
coordinate  system  and  then  map  the  results  to  the  latitude-longitude  system.  The  definition  of  the 
auxiliary  coordinates,  which  are  denoted  here  by  (0,y),  is  illustrated  in  Fig.  2.  In  this  definition,  it  is 
assumed  that  the  route  set  is  such  that  there  is  a  “nominal  route”  co„  with  the  property  that,  for  each  point 
(A/,  4/)  on  co„ ,  there  is  a  great  circle  arc  through  that  point  that  is  orthogonal  to  the  nominal  route  and 
that  intersects  every  other  route  in  the  route  set.  The  auxiliary  coordinate  0  is  defined  as  the  distance 
along  the  nominal  route  from  the  departure  port  to  the  point  (A/,  4/) .  This  distance  increases  from  zero  at 
the  departure  port  to  Ln  at  the  destination  port,  where  Ln  is  the  length  of  the  nominal  route  expressed  in 
radians.  The  “cross-sectional  coordinate”  y  is  defined  as  the  signed  distance  along  the  great  circle  arc 
from  the  point  (A/,  4/)  on  the  nominal  route  eo„  to  the  point  ( /. " ,  cj) " ) .  This  distance  is  taken  to  be  positive 
in  the  direction  of  the  upper  route  envelope,  denoted  by  eu(Q )  in  Fig.  2,  and  negative  in  the  direction  of 
the  lower  route  envelope  e/(9).  For  route  segments  that  lie  in  the  virtual  coordinate  set, 
(A,  4>)  e  sjf  v  =  (-go, — tc)  x  £2 ,  we  take  0  =  A  +  n  and  y  equal  to  a  function  of  <|)  that  is  specified  below. 


Ship-Track  Models  Based  on  Poisson-Distributed  Port-Departure  Times 


5 


Fig.  2  —  The  auxiliary  coordinate  system 


With  the  auxiliary  coordinate  system  so  defined,  the  longitude-latitude  coordinates  are  related  to  the 
auxiliary  coordinates  by  a  one-to-one  transformation.  Consequently,  the  route  function  Rf  ^  (x;  eo)  can  be 
obtained  by  first  determining  this  function  in  terms  of  the  auxiliary  coordinates  (i.e.,  as  the  function 
Rq  y  (x;  eo) )  and  then  using  the  coordinate  transformation  to  express  the  results  in  the  longitude-latitude 
coordinate  system.  Reference  5  provides  an  example  of  a  one-to-one  transformation  between  the  auxiliary 
coordinates  and  the  longitude-latitude  coordinates  for  a  particular  choice  of  the  nominal  route. 


In  the  deterministic  model  presented  in  Section  4,  all  of  the  routes  in  the  route  set  “progress”  from  the 
departure  port  to  the  destination  port  without  “doubling  back”  on  themselves.  For  this  case,  there  can  be 
only  one  value  of  the  cross-sectional  coordinate  y  for  each  value  of  the  nominal  route  coordinate  9 . 
Consequently,  for  each  route  eo,  the  cross-sectional  variable  can  be  expressed  as  a  function  y  =  a(9,eo)  , 
so  that  the  route  eo  is  traced  out  as  the  points  (9,a(9,co))  as  9  increases  from  zero  to  Ln  .  Furthermore, 
the  distance  traveled  along  the  route  is  given  by 


X(9;eo) 


1  + 


5a(9',oo) 


X2 


-il/2 


59' 


d  9' 


(3) 


which  is  an  increasing  function  of  the  nominal  route  coordinate  9 .  ft  follows  that  9  can  be  expressed  as  a 
function  of  the  distance  traveled  by  the  inverse  function 

0(x;eo)  =  X~l  (x;eo)  .  (4) 

With  these  definitions,  the  route  function  for  progressive  route  sets  is  given  by 

Rq  (x;co)  =  0(x;eo)  =  X~l  (x;eo) 

Ry  (x;co)  =  a(0(x;eo),eo) 


(5) 
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2.2  The  Ship-Motion  Function 

The  ship-motion  function  describes  the  distance  traveled  by  a  ship  along  a  route  as  a  function  of  the 
elapsed  time  after  its  departure  from  the  port.  Since  this  motion  depends  on  the  nominal  speed  that  the 
ship  travels  and  possibly  on  the  specific  route  that  the  ship  follows,  we  include  these  quantities  as 
parameters  of  the  ship-motion  function.  Furthermore,  we  also  allow  the  ship-motion  function  to  depend 
on  absolute  time  to  account  for  environmental  effects  (e.g.,  storms).  With  these  conventions,  the  distance 
traveled  from  the  departure  port  at  time  t  by  a  ship  that  departs  at  time  t  and  then  travels  the  route  co 
with  nominal  speed  v  is  represented  by  x  -  M(t  -  x;f,co,v) .  As  defined,  the  ship-motion  function  takes 
on  distance  values  in  the  interval  [0,  L(co)]  as  t  varies  from  the  departure  time  x  to  the  arrival  time  at  the 
destination  port  rj .  As  with  the  route  function,  it  is  convenient  to  allow  the  distance  traveled  to  take  on 
“virtual”  values  by  extending  the  definition  of  the  ship-motion  function  to  the  whole  time  axis.  To  this 
end,  we  assume  constant  speed  motion  for  t  £  [x,r|]  and  define  an  extended  ship-motion  function  by 

v(f-x)  fort<x 

M(V-x;t,co,v)  forte[x,r|]  .  (6) 

v(t-x)  +  L(co)  fort>x 

With  this  definition,  the  extended  motion  function  Mx(t-x;t,  co,v)  is  a  monotonically  increasing 
function  of  elapsed  time  q  =  t  -  x  .  It  follows  that  given  any  initial  time  t0 ,  there  is  a  unique  initial 
distance  x0  obtained  by  setting  q  =  q0  =  t0  -  x  in  Eq.  (5).  Furthermore,  given  an  initial  distance  x0  ,  there 
is  a  unique  value  of  the  initial  elapsed  time 

*30  ( X0  ’  Aj  )  —  t0  _  t  —  M X  (X0  To  5  ®>v)  i  (7) 

corresponding  to  x0  .  Consequently,  it  is  always  possible  to  express  the  distance  traveled  as  a  function  of 
the  initial  time  t0  and  the  initial  distance  x0  .  To  this  end,  we  redefine  the  ship-motion  function  by 

M  (t -t0;x0,(o,v)  =  Mx  (t-ta  +q0  (x0,co,v;t0);t,(o,v)  .  (8) 


Mx(t-v,t,  oo,v)  =  j 


Figure  3  shows  an  example  of  these  definitions.  In  this  figure,  we  have  illustrated  the  ship-motion 
function  for  two  ships,  each  of  which  travels  the  same  route  co  with  the  same  nominal  velocity  v .  The 
first  ship  departs  at  time  Xj  and  arrives  at  time  //, ;  the  second  departs  at  time  x2  and  arrives  at  time  q2  • 
During  the  period  when  the  second  ship  is  in  transit,  a  local  storm  causes  it  to  reduce  speed  with  the  result 
that  its  total  transit  time  q2  -  x2  is  longer  than  the  transit  time  for  the  first  ship.  Both  ships  are  present 
during  the  observation  interval  [t0,t0+T].  The  first  ship  is  en  route  at  the  initial  time  t0  since 
t0  e(xj,r|1).  This  ship  has  an  elapsed  time  q(x0 ,  t0 ; co,  v)  >  0  and  a  “real”  initial  distance  x0  e[0,T(co)]. 
The  second  ship  in  not  en  route  at  the  initial  time  since  t0  e  (x2,q2) .  However,  it  is  en  route  during  the 
observation  interval  since  t0  +T  e  (x2,q2) .  This  ship  has  an  elapsed  time  q(x0,t0; co,v)<0  and  a 
negative  initial  distance  x0  e  [vT,  0] .  The  first  ship  arrives  at  the  destination  port  in  the  observation 
interval;  the  second  ship  does  not. 
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Note  that  in  this  example,  the  nominal  speed  parameter  specifies  a  “characteristic”  speed  of  the  ship 
rather  than  the  actual  speed  made  good.  Other  attributes  of  the  ship  relevant  to  a  particular  application  of 
the  model  (e.g.,  type  of  ship)  can  be  included  in  the  model  by  interpreting  v  as  a  vector  parameter  in  the 
subsequent  development.  Also  note  that  the  definition  of  the  motion  function  does  not  preclude  the 
second  ship  from  overtaking  and  passing  the  first  ship.  Strictly  speaking,  with  the  route  function  defined 
by  Eq.  (1),  this  would  amount  to  the  second  ship  “passing  through”  the  first  ship  if  both  ships  were 
traveling  on  the  same  route.  However,  with  the  specific  route  functions  considered  in  this  report,  this  is  an 
event  of  zero  probability  and  is  not  considered  further. 

2.3  The  Ship-Track  Function 

The  ship-track  function  Gy  ^  is  obtained  from  the  route  function  Ry  ^  by  using  the  ship-motion 
function  M  to  express  the  distance  traveled  by  a  ship  as  a  function  of  the  elapsed  time  and  the  ship 
parameter.  It  follows  from  Eqs.  (1)  and  (7)  that  the  longitude  coordinate  function  Gy  and  the  latitude 
coordinate  function  G§  of  the  ship-track  function  have  the  form 

Gx(t  -t0;x0,(o,v)  =  RX(M (t -t0;x0,(o,v); co) 

Gty  ( t  - 10 ;  x0 ,  co,  v)  =  [M  ( t  -  tQ ;  x0 ,  co,  v) ;  co) 

Note  that  only  the  route  parameter  component  of  the  track  parameter  ( x0 ,  co,  v)  is  needed  to  identify  the 
particular  route,  whereas  the  complete  track  parameter  is  needed  to  specify  the  motion  along  that  route. 

A  ship-track  realization  is  obtained  by  specifying  the  number  of  ship  tracks  in  the  realization  n  and 
the  track  parameters  for  those  ships  {(x0k  ,co^,v^);  k  =  1 ,...,«}  and  then  computing  the  track  of  each  ship 
using  the  track  function.  An  example  of  a  ship  track  realization  is  shown  in  Fig.  4.  The  routes  that  the 
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ships  travel  are  shown  by  the  thin  curves.  The  tracks  of  the  ships  along  those  routes  for  the  time  interval 
[tQ,t0  +  7’]  are  indicated  by  the  heavy  lines  on  these  curves.  For  the  tracks  along  the  routes  coj  and  CO4, 
the  ships  are  en  route  during  the  whole  time  interval  \t0,t0  +  7  ] .  For  the  track  along  the  route  0)3,  the 
ship  is  en  route  at  time  t0  but  it  arrives  at  the  destination  port  prior  to  ta+T  .  Thus,  for  these  three  tracks, 
the  initial- distances  are  positive  and  the  initial  positions,  indicated  by  large  dots,  lie  in  the  set  of  real  ship 
coordinates.  For  the  track  along  the  route  co2  >  the  ship  departs  from  Pn  after  time  t0  .  Therefore,  the  initial 
distance  for  this  track  x0l  is  negative  and  the  initial  position  lies  in  the  set  of  virtual  ship  coordinates.  At 
time  t0 ,  the  ship  starts  at  the  initial  distance  x0n  and  moves  along  the  route  segment  in  the  virtual  ship 
coordinates  with  constant  speed  V2  up  to  the  departure  time  x2  =t0  -  x(n  /v2  .  At  this  time,  it  actually 
departs  the  port  and  travels  along  the  route  segment  shown  with  its  motion  determined  by 
M(t-t0;t,x02,  co2,v2). 


Fig.  4  —  Example  of  a  ship  track  realization  for  a  single  port-pair 

To  determine  individual  track  realizations,  it  is  necessary  to  have  a  mechanism  for  determining  the 
track  parameters  for  each  realization.  This  is  done  in  terms  of  the  port  parameter  probability  law  defined 
in  the  following  section. 

3,  THE  PORT-DEPARTURE  AND  TRACK-PARAMETER  PROBABILITY  LAWS 

To  an  observer  located  at  the  departure  port,  the  shipping  is  described  by  the  sequence  of  parameters 
{(t£,oc>£,V£)}  ,  where  is  the  departure  time  of  the  kth  ship  to  leave  the  port,  X/.  is  its  route  parameter, 
and  is  its  nominal  speed.  In  this  report,  we  assume  that:  ( 1 )  the  ship  departure  times  are  described  by  a 
Poisson  process  with  a  time-dependent  rate  function  pT  (?)  that  represents  the  mean  number  of  ship 
departures  per  unit  time;  (2)  the  route  that  each  ship  travels  and  its  speed  are  independent  of  its  departure 
time  and  independent  of  one  another;  (3)  different  ships  select  routes  independently  of  one  another  and 
each  ship  selects  its  route  from  the  same  probability  density  pM  ( co) ;  and  (4)  the  nominal  speeds  of  all 
ships  are  independent  of  one  another  and  are  described  by  the  same  probability  density  pv  (v) .  The  rate 
function  pT(t)  can  be  estimated  from  ship  departure  data;  the  ship  speed  density  pv  (v)  can  be 
determined  from  the  ships’  register  data  given  the  identities  of  the  departing  ships.  However,  in  general, 
the  route  parameter  probability  density  pm(co)  cannot  be  determined  independently  of  the  probability  law 
describing  the  route  set. 

According  to  the  Poisson  departure  time  assumption,  the  number  of  ships  leaving  the  port  during  the 
time  interval  (ta ,to+T]  is  a  Poisson  random  variable  Nx  with  a  probability  mass  function  given  by 
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Pr  [Nt 


=  „]  =  exp{-Mt((0,r)}^%3 

n ! 


(10) 


where 


Mx(t0,T)  =  l‘0+TpJt)dt 

lO 


(11) 


is  the  mean  number  of  ships  to  leave  the  port  in  the  interval.  Furthermore,  given  that  there  are  n  ship 
departures  during  the  interval,  the  ship  departure  times  \xk  ;k-\,...,n}  are  independent  and  identically 
distributed  with  common  probability  density. 


pT(xkU0J) 


V-x(ik)Mx{to>T)  1  forTt  e[t0’to+T] 
0  otherwise 


(12) 


Note  that  since  the  departure  rate  px(0  depends  on  time,  the  mean  number  of  ships  Mx  (t0 ,  T) ,  the 
departure-time  probability  density  px (xk  ;t0,T) ,  and  the  probability  mass  function  depend  on  the  choice 
of  the  time  interval  \t0,t0  +  7] .  We  allow  for  this  time  dependence  to  incorporate  variations  in  the  number 
of  ship  departures  with  time  of  day,  day  of  the  week,  etc.  For  homogeneous  Poisson  processes,  the 
departure  rate  is  independent  of  time  (pT(t)  =  pT);  hence,  the  mean  number  of  ship  departures  is  simply 
T pT  and  the  departure  times  are  uniformly  distributed  on  the  interval  \tQ ,ta+T]  (i.e.,  pT(xk)  =  T_1 ). 

It  follows  from  the  port-departure  assumptions  that  the  port-departure  parameter  (x,co,v)  is  described 
by  a  multidimensional  Poisson  process  that  is  specified  by  a  rate  function  pT  a  V(t,oo,v)  .  Moreover,  from 
the  definitions  of  Section  2,  there  is  a  one-to-one  transformation  between  the  port-departure  parameter 
(i,co,v)  and  the  ship-track  parameter  (x0,co,v) .  It  follows  that  the  ship-track  parameter  is  also  distributed 
as  a  multi-dimensional  Poisson  process  that  is  specified  by  a  rate  function  px  mv(x0,(o,v;t0,T) .  The 
specification  and  derivation  of  these  rate  functions  is  presented  in  Appendix  B  (see  Eqs.  (B 1),  (B4),  and 
(B5)).  The  definition  of  a  multidimensional  Poisson  process  is  presented  in  Appendix  A. 

The  Poisson  process  on  the  ship-track  parameters  describes  the  number  and  the  distribution  of  the 
parameters  for  all  ships  present  in  the  region  during  the  interval  [t0 ,t0  +T],  As  seen  in  Section  2,  the 
ships  present  in  the  region  during  this  interval  are  those  that  are  en  route  at  time  t0  and  those  that  depart 
during  the  half  open  interval  (t0,t0+T\.  The  former  have  positive  initial  distances;  the  latter  have 
negative  initial  distances.  Consequently,  by  restricting  the  track-parameter  rate  function 
Px  ,o),v(xo’ wAT0,T)  to  positive  values  of  x0  ,  we  obtain  a  Poisson  process  that  describes  the  ships  that 
are  en  route  at  time  t0 .  Similarly,  by  restricting  PXq ,co  v(xo,0),v;to,T)  to  negative  values  of  x0  ,  we  obtain 
a  process  describing  the  ships  that  depart  during  (t0,t0  +  7  ] .  The  total  number  of  ships  en  route  during 
[to,t0  +T]  is  the  Poisson  process  obtained  as  the  composition  of  the  two  component  processes.  For  each 
of  these  processes,  the  number  of  ships  is  a  Poisson  random  variable  with  a  probability  mass  function 
given  by  Eq.  (10)  and  the  track  parameter  for  those  ship  are  independent  and  identically  distributed  with  a 
probability  density  given  by  an  equation  analogous  to  Eq.  (12).  The  mean  number  of  ships  is  obtained  as 
an  integral  of  pYo  v  (x0,eo,v;70,r)  over  the  appropriate  volume.  The  definition  of  these  processes  and 
explicit  formulas  for  the  relevant  quantities  are  presented  in  the  following  sections. 
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The  track-parameter  rate  function  px  (av(x0,a> ,v;t0,T)  completes  the  model  if  the  route  set 
probability  density  pa( eo)  and  the  route  function  of  Eq.  (1)  are  known.  In  a  simple  approximate  model, 
where  the  traffic  is  described  in  terms  of,  at  most,  a  finite  number  of  routes,  the  route  function  might  be 
described  analytically  in  terms  of  great  circle  route  segments  and  the  weightings  pa( co)  might  also  be 
available.  However,  in  the  more  realistic  case,  where  there  is  a  continuum  of  routes,  the  route  function 
and  the  route  parameter  density  must  be  determined  from  data  on  the  routes  that  the  ships  follow.  In 
principle,  these  data  can  be  used  to  construct  a  probability  law  on  the  ensemble  of  routes  from  Pn  and 
Pm  .  In  general,  such  a  probability  law  consists  of  the  joint  probability  density  on  the  coordinate  positions 
{(}i(x1),())(xi)),...,(^(x„),(t)(x„))}  for  all  possible  samples  of  the  distances  traveled  {jq  <  ...<  xn;  n  >  0} . 
To  determine  such  a  probability  law  requires  more  data  than  are  usually  available.  Consequently,  it  is 
necessary  to  introduce  approximations  and  assumptions  that  allow  the  route  functions  to  be  defined  in 
terms  of  simplified  probabilistic  descriptions. 

In  the  following  section,  we  describe  a  “deterministic”  model  where  the  route  functions  are 
determined  from  a  probability  density  that  can  be  estimated  from  any  sample  of  ship  coordinates, 
regardless  of  whether  those  coordinates  are  organized  into  specific  routes.  In  Section  5  we  describe  a 
“stochastic”  model  which  assumes  that  the  ship  routes  are  described  by  a  Markov  process. 

4.  A  DETERMINISTIC  ROUTE  FUNCTION  MODEL 

In  the  ship-track  model  presented  here,  both  the  route  function  and  the  route  parameter  density  are 
determined  from  a  probability  density  on  the  cross-sectional  ship  positions  y  conditioned  on  the  nominal 
route  coordinate  positions  0 .  This  is  done  in  such  a  way  that  each  route  is  uniquely  determined  by  any 
point  on  that  route.  As  a  consequence,  the  Poisson  process  describing  the  distribution  of  ship-track 
parameters  can  be  expressed  in  terms  of  the  initial  coordinates  of  the  ship  at  time  t0 ,  rather  than  the  initial 
distance  and  the  route  parameter.  This  simplicity  is  obtained  by  limiting  the  route  structure  to  progressive 
route  sets  with  the  additional  constraint  that  the  individual  routes  can  not  intersect  one  another.  To 
describe  the  model,  we  first  define  the  underlying  probability  densities.  We  then  specify  the  route  and  the 
track  functions  that  are  determined  by  these  densities.  Finally,  we  present  the  probability  law  on  the  track 
parameter  which  results  from  these  definitions  and  the  assumptions  of  Section  3. 

4.1  The  Route  Coordinate  Probability  Density 

The  route  coordinate  probability  density  describes  the  concentration  of  the  y  values  of  the  routes  as 
they  intersect  the  cut  in  the  route  set  at  0.  We  denote  this  probability  density  by  /i.,q  (y;  0 J ,  where  0 
specifies  the  location  of  the  cut,  and  refer  to  the  corresponding  cumulative  probability  distribution 

y 

f  ^yie(Y';e)^y'  (13) 

e/(0) 


as  the  route  coordinate  distribution  function.  Figure  5  shows  an  example  of  a  route  coordinate  density  and 
distribution  function  for  a  specific  cut  in  the  route  set.  In  this  example,  the  routes  for  the  cut  0  =  0O ,  are 
concentrated  near  the  lower  envelope  e/(0o).  The  route  coordinate  density  at  other  cuts  in  the  route  set 
need  not  be  the  same  as  the  one  shown.  For  example,  the  distribution  of  routes  shown  in  the  figure 
suggests  that  the  route  coordinate  density  for  the  cut  at  0  =  0j  would  reflect  a  more  disburse 
concentration  of  the  routes  closer  to  the  center  of  the  route  envelopes.  A  method  for  estimating  these 
densities  from  ship  route  data  can  be  found  in  Ref.  2  along  with  a  number  of  examples. 
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upper  route  route,  m 


Fig.  5  —  Example  of  a  route  coordinate  probability  density 

It  is  important  to  emphasize  that  the  route  coordinate  densities,  in  themselves,  do  not  specify  the 
routes  that  the  ships  travel  since  many  different  possible  route  functions  can  satisfy  a  given  route 
coordinate  density. 

4.2  The  Route  Function 

For  progressive  routes,  the  route  sets  are  specified  by  the  cross-sectional  function  a(0,co)  and  the 
route  parameter  probability  density  pbJ(co) .  For  the  deterministic  model,  we  define  these  quantities  in  such 
a  way  that  the  resulting  set  of  routes  is  consistent  with  the  route  coordinate  distribution  function  (y;  9)  . 
To  this  end,  we  restrict  eo  to  the  interval  [0,1]  and  define  the  cross-sectional  function  a(0,ro)  as  the 
inverse  of  the  route  coordinate  distribution.  Specifically,  for  each  eo  e  [0,1] ,  a(0,eo)  is  the  function 
satisfying 

Py|6(a(0,co);0)  =  co  .  (14) 

Clearly,  a(9,ro)  is  well  defined  since  for  each  value  of  0,  (y;  0)  is  an  increasing  function  of  y 

taking  values  in  the  interval  [0,1];  hence,  for  fixed  0,  there  is  an  inverse  function  that  maps 
roe[0,l]  to  the  set  of  all  possible  cross-sectional  coordinate  values  y.  For  a  fixed®,  the  y  value 
corresponding  to  the  0  value  on  the  route  ro  is  the  value  given  by  this  inverse  function.  For  notational 
convenience,  we  write 

a(0,ro)  =  (ro;0)  for  0>0  .  (15) 

Figure  6  illustrates  the  definition  of  the  cross-sectional  function.  As  seen  in  the  figure,  the  two  route 
envelopes  are  also  routes  in  the  route  set.  The  lower  route  envelope  e/(0)  corresponds  to  the  route 
parameter ®  =  0 ;  the  upper  route  envelope  eu(Q)  corresponds  toro  =  l.  For  fixed  0,  the  cross-sectional 
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function  a(0,co)  yields  y  values  that  increase  over  the  interval  [cy  (0),et((0)]  as  co  increases  over  the 
interval  [ 0,  l] .  For  a  fixed  value  ofeo,  a(0,co)  traces  out  the  y  values  along  the  route  co  as  0  increase 
from  0  to  Ln  .  Note  that  the  routes  are  defined  so  that  no  two  can  intersect  one  another.  (If  they  could,  Eq. 
(13)  would  not  be  well-defined.)  Also  note  that  since  the  lower  route  envelope  and  the  upper  route 
envelope  are  assumed  to  be  distinct,  the  interval  [e/(0),eM(0)]  cannot  degenerate  to  a  single  point  even  at 
the  departure  port  0  =  0  or  at  the  destination  port  0  =  L  . 


Fig.  6  —  Definition  of  the  cross-sectional  route  function  and  the  route  parameter 


Equation  (15)  describes  the  segments  of  the  routes  that  lie  between  the  departure  port  and  the 
destination  port.  These  route  segments  lie  in  the  set  of  real  ship  coordinates  SJ1;.  =  [O,Z.]x[^(0),eu(0)]. 
For  the  route  segments  that  describe  the  fictitious  motion  of  the  ships  before  they  leave  the  port,  we  take 
a(0,co)  to  be  the  inverse  route  distribution  function  evaluated  at  the  departure  port  (i.e., 
a^co^P^tco^O))  for  0  <  0 .  These  route  segments  lie  in  the  set  of  virtual  ship  coordinates 
iHv  =  (-co,  0)  x  [cy  (0),  eu  (0)]  and  have  y  coordinates  that  are  independent  of  0 . 

The  definition  of  the  cross-sectional  function  suffices  to  specify  the  route  function  Rq  (x;  co)  .  The  0 
coordinate  of  Rq^(x;  co)  is  @(x0;co)  =  X_1(x;co),  where  the  distance  traveled  function  X(0;co)  is 
determined  from  a(0,co)  by  Eq.  (3).  The  y  coordinate  of  Rq^(x;  co)  is  a(@(x;co),co) .  The  track  function 
Gq  y(t  -t0;x0,(o,v)  is  obtained  from  the  route  function  using  the  motion  function  M(t-t0;x0, co,v)  to 
express  the  distance  traveled  as  a  function  of  time. 

An  important  consequence  of  the  definition  of  the  cross-sectional  function  is  that  there  is  a  one-to-one 
transformation  between  the  track  parameter  components  (x0 ,  co)  and  the  initial  coordinates  (0O ,  y0 ) .  In 
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particular,  the  initial  coordinates  are  determined  by  90  =  0(xo;co)  andy0  =  a(0(xo;co);co) .  The  initial 
distance  and  the  route  parameter  are  determined  by 


«=^o(Yo;O)=^y!e(Yo;0) 


(16) 


and 


'■(J 

xo  ~  *0  (0o’Yo)  ~  J" 


-.1/2 


1  + 


dPy\Q~]  (“;0') 


50' 


c/9'  . 


(17) 


Equation  (16)  follows  immediately  from  Eq.  (14)  by  replacing  a(9,eo)  by  y0  .  Equation  (17)  follows 
from  Eq.  (3)  by  setting  9  =  90  and  by  substituting  for  eo  from  Eq.  (15).  Note  that  by  virtue  of  Eqs.  (16) 
and  ( 1 7),  the  route  function  can  be  expressed  in  terms  of  the  initial  position  and  the  motion  function  and 
the  track  function  can  be  expressed  in  terms  of  the  initial  position  and  the  characteristic  speed. 
Consequently,  the  ship  track  realizations  can  be  obtained  from  the  track  function  by  specifying  the  track 
parameter  sets  {(0O  ,y0  ,  vy. );  k  - 1 ,...,«} .  The  interpretation  of  a  ship  track  realization  is  the  same  as  that 

of  Fig.  3  with  x0  and  co  replaced  by  90  and  y0  (i.e.,  the  initial  positions  (the  dots)  are  the  points 
(90,y0)  and  the  route  parameter  in  the  ship-motion  function  is  determined  from  these  points  by  Eq. 
(16)).  The  probability  law  on  the  track  parameters  is  described  in  the  following  subsection. 

To  conclude  these  definitions,  we  show  that  the  cross-sectional  function  is  consistent  with  the  route 
coordinate  distribution  function  if  the  route  parameter  is  uniformly  distributed.  To  see  this,  we  note 
that  ^,|e  (y;  9)  =  Prob{y'  e  [e/(9),y]}  =  Prob{y'  e  [ey  (9),a(9;eo)]} .  But  y'  e  [e/(9),a(9;eo)]  is  equivalent 

toco'e[9,Ey|0(y;9)],  so  thatProb{y' e [e?/ (0), a(9; oo)] }  =  Prob{coe [O,Py|0(y;0)]} ,  which  is  the  same  as  the 

cumulative  distribution  of  the  route  parameter  evaluated  at  eo  =  Ey|0(y;9) .  Thus,  E,|0(y;O)  =  Eco(/^0(y;0)) , 

where  P(l)  ( co)  is  the  route  parameter  cumulative  distribution  function  determined  from  the  route 
parameter  probability  density  pa  ( co) .  This  equality  holds  if  and  only  if  the  route  parameter  is  uniformly 
distributed  on  [0, l] ,  in  which  case  P(,, ( co)  =  co . 

4.3  The  Track-Parameter  Process 

The  ship-track  probability  law  describes  the  number  and  the  distribution  of  the  parameters  of  all  the 
ships  present  in  the  region  during  the  interval  \t0,t0  +  7] .  For  the  deterministic  model,  this  probability 
law  is  obtained  from  the  probability  law  of  Section  3  through  an  invertible  transformation  that  maps  the 
parameter  ( x0 ,  co,  v)  to  the  parameter  (90 ,  y0 ,  v) .  It  is  shown  in  Appendix  B  that  this  results  in  a  Poisson 
process  with  a  rate  function  p0,y,v(Qo,y 0,v;l0,T)  defined  on  the  extended  route  set  utR,. .  As 

noted  in  Section  3,  this  process  can  be  represented  as  the  composition  of  two  Poisson  processes,  one  that 
describes  the  ships  that  are  en  route  at  time  t0  and  one  that  describes  the  ships  that  depart  during 
(t0,t0  +  T]  ■  We  refer  to  the  former  as  the  “en  route”  process  and  to  the  later  as  the  “entry”  process. 

To  define  these  processes,  it  suffices  to  specify  their  rate  functions.  The  number  of  ships  and  the 
probability  density  on  the  parameters  for  those  ships  are  then  given  by  Eqs.  (10)  and  (12),  where  the  mean 
number  of  ships  is  obtained  as  an  integral  over  the  rate  function.  To  compute  this  integral,  it  is  convenient 
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to  define  a  shipping  density  Pq  y  (90,y0;t0)  as  the  integral  of  the  rate  function  with  respect  to  ship  speed. 
The  mean  number  of  ships  is  then  obtained  as  the  integral  of  the  shipping  density  with  respect  to  the  ship 
coordinates.  The  shipping  distribution  is  the  integral  of  the  shipping  density  over  element  areas.  The 
specific  equations  are  as  follows.  The  derivations  of  these  equations  are  presented  in  Appendix  B. 

For  the  en  route  process,  the  only  ships  that  are  en  route  at  time  t0  are  those  that  have  real  ship 
coordinates.  Consequently,  the  rate  function  is  obtained  by  restricting  pQ y  V(90,y0,v;t0,r)  to  the  real 
ship  coordinates  iH,.  =  [9,Z]x[e;(9),e„(9)] .  The  rate  function  is  given  by 


Fo,y,v  (®o> Yo  ATo  )  ■ 


Fx  {fo  )) 


(  x0 , 0),  v,  t0  ) 


dxn 


Xo=Xo(Qodo) 

®=^o(Yo’00) 


Py\Q  (Vo  ) 


s^0(e0,Y0) 


59, 


Pv{v) 


(18) 


As  noted  above,  the  shipping  density  for  the  ships  that  are  en  route  at  time  t0 ,  Pq  ., (90 , y 0;t0)  is  the 
integral  of  the  rate  function  with  respect  to  ship  speed; 

00 

M-9,y  (00’  Yo To  )  —  j"  FO,y,v  (®o  > Yo  ATo ) dv,  (90,y0 )  £  .  (19) 

0 

The  mean  number  of  ships  with  initial  coordinates  (A.',  ())r)  in  any  set  Cr  c  91,.  is  the  integral  of  the 
shipping  density  with  respect  to  the  ship  coordinates  over  the  set  Cr ; 

Mr{t0\Cr)  =  jc  P0,y (9',y';t0)t/9'i/y'  .  (29) 

It  follows  from  the  definition  of  a  Poisson  process,  that  the  number  of  ships  en  route  at  time  t0  with 
coordinates  (90,y0)  e  Cr  c  $Rr ,  Nr(t0;Cr),  is  a  Poisson  random  variable  with  a  probability  mass 
function  given  by  Eq.  (19)  with  M x  replaced  by  Mr(t0;Cr) .  Furthermore,  the  track  parameters  for  those 
ships  are  independent  and  identically  distributed  with  the  common  probability  density  function 

/V (Ho ’Yo aTq > C;. )  — p0y>v (90 ,y0 ,v,f0 (Y0 , C,. )  (90,y0)GtH?.  .  (21) 

The  total  number  of  ships  en  route  at  time  t0  is  determined  from  these  equations  by  taking  Cr  -  . 

An  interpretation  of  these  equations  is  as  follows.  The  first  factor  in  Eq.  (18)  describes  the  mean 
number  of  ships  per-unit  distance  along  the  route  determined  by(90,y0) .  This  factor  incoiporates  the 
temporal  variations  in  the  rate  at  which  ships  depart  the  port  through  the  time  dependence  on  the  ship 
departure  rate  pT(t)  and  the  motion  of  the  ships  along  the  route  through  its  dependence  on  the  departure 
time  function  go  (xo,  a>,  v;  ta ) ,  and  its  first  derivative.  For  the  special  case  of  constant  speed  motion, 
go(xo,C0,v;to)  =  xov~l ,  this  factor  simplifies  to  pT(t0  -  X0(90,y0)v-1)v-1 .  Thus,  for  constant  speed 
motion,  a  change  in  the  rate  at  which  ships  depart  the  port  simply  propagates  along  the  route  with  the 
nominal  ship  speed. 

The  second  factor  in  Eq.  (18)  describes  how  the  ships  are  distributed  across  the  routes.  This  factor  is 
equal  to  the  route  coordinate  density  weighted  by  the  magnitude  of  a  derivative  that  represents  the  rate  of 
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change  in  the  distance  traveled  with  respect  to  the  nominal  route  distance.  For  most  route  sets,  the  routes 
do  not  depart  significantly  from  the  nominal  route  so  that  this  derivative  is  approximately  unity  and 
hence,  the  second  factor  is  approximately  equal  to  the  route  coordinate  density.  For  these  “narrow”  route 
sets  and  for  constant  speed  motion,  the  track-parameter  rate  function  simplifies  to 

F0,y,v  (®o  )  —  Ft  i^o  ~  (0o  >Yo  ("/o’^o  )1;  Pviy)>  (0<3’Yo  )  e  ’  (22) 

and  the  shipping  density  becomes 

00 

F0,y  (®o  ’  y  o’^o )  —  Fyj0  (Yo  )  J"  Ft  i^o  ~~  (0o ,  Yo  )v  ji;  Pv{y)dv,  (^0o,yo)gJ?;.  .  (23) 

0 

Note  that  if  in  addition  the  departure  rate  is  independent  of  time,  pT(x)  =  pT,  then  the  rate 
function  is  p0y>v(0o,Yo^y;to,r)  =  pTpy|0(Yo;0o)v_1/>v(v)  and  the  shipping  density  is  p0y (0o,Yo;to)  = 
FxPPy|e(Yo;0o),  where  P  is  the  mean  value  of  the  reciprocal  speed,  $  =  E[v  *]. 

For  the  entry  process,  the  only  ships  that  depart  in  the  interval  (to  ,ta+T]  are  those  that  have  virtual 
ship  coordinates.  Consequently,  the  rate  function  for  the  entry  process  is  obtained  by  restricting 
p0  y  v(0o,Yo,v;to,r)  to  the  virtual  ship  coordinates  iRv(r)  =[-vF,0]x  [^(0  =  O),eJ((0  =  0)] .  The  rate 
function  is  given  by 

lle,y;v(0o>YO)v;fo,7’)  =  pT(to-0ov"1)py|0(Yo;0o  =0)v_1j9v(v);  (0o,yo)  GiRv(r),  (24) 
and  the  shipping  density  is 

00 

F0,y(0o>Yo;^>r)  =  Py|e(Yo;0o  =0)  j  \ix(t0-Q0v~iy~lpv{v)dv,  (0o,Yo)G0?v(r)  .  (25) 

K\/t 

Note  that  Eqs.  (24)  and  (25)  are  the  same  as  Eqs.  (22)  and  (23),  except  that  _py|0(Yo;9o)  is  replaced 
by  py|9(Y0;90  =  0)  and  the  lower  limit  of  the  integral  in  Eq.  (24)  depends  on  both  (/.',<)')  and  the  time 
interval  duration  T .  This  is  not  surprising  since,  by  the  definition  of  the  route  function  on  the  virtual 
coordinates,  the  ships  move  at  constant  speed  and  the  derivative  of  XQ  ( 0o  ;y0)  is  unity.  Also  note  that  for 
the  entry  process,  the  rate  function  and  the  shipping  density  depend  on  the  duration  T  of  the  interval 
\to,t0  +  T] ,  as  well  as  the  initial  time  t0 ;  whereas,  for  the  en  route  process,  these  parameters  depend  only 
on  t0. 

Equations  (24)  and  (25)  determine  the  entry  process.  The  mean  number  of  the  ships  that  depart  in  the 
interval  (t0,t0  +T]  with  coordinates  (0o ,yr) )  e  Cv  c  '.)?,,( F ) ,  M T(t0,T; Cv ) ,  is  obtained  by  integrating  the 
shipping  density  over  the  set  Cv  ctRv(F) .  The  probability  mass  function  on  the  number  of  these  ships 
Ne(t0,T;Cv )  is  given  by  Eq.  (9)  with  MT  replaced  byMe(t0,F;Cv).  The  probability  density  on  the 
parameters  of  these  ships  pe(Q0,y0,v,t0,T)  is  given  by  Eq.  (11)  with  pT(x^)  replaced  by 
F0,y,v(0O’YoA;?o,F)  and  Mx(t0J)  replaced  by  Me(t0,T :  Cv)  . 

Finally,  the  Poisson  process  describing  the  total  number  of  ships  in  the  interval  Np(t0,T )  is  the 
composition  of  the  en  route  and  the  entry  processes.  Thus,  for  the  closed  interval  \t0 ,  t0  +  T\ ,  the  total 
number  of  ships  is  NT  (t0 ,  T)  =  Nr  (t0 )  +  Ne  (t0 ,  T) ,  the  mean  number  of  ships  is 
Mt  (t0 ,  T)  =  Mr  (t0 )  +  Me  (t0 ,  T) ,  and  the  track-parameter  probability  density  is 

PT (0O . Y0 , v; t0 , T)  =  n0>y >v (0O , Y0 ,v, t0 , T) / MT (tQ , T)  . 
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The  route  set  assumptions  for  the  deterministic  model  impose  potentially  important  restrictions  on  the 
routes  in  the  route  set.  In  particular,  the  route  set  approximation  can  hold  only  if  the  routes  have  no 
common  points  and  only  if  individual  routes  do  not  have  multiple  points  (i.e.,  the  routes  cannot  intersect 
one  another  and  cannot  cross  over  themselves).  For  those  ships  that  simply  transit  from  one  port  to 
another,  the  approximation  that  results  from  these  restrictions  may  be  adequate  for  those  applications  that 
are  limited  to  port-to-port  traffic  and  do  not  require  precise  descriptions  of  the  tracks.  There  is,  however, 
shipping  traffic  which  does  more  than  simply  transit  from  one  port  to  another  (e.g.,  fishing  vessels, 
recreational,  and  military).  For  applications  where  this  traffic  is  important,  the  route  set  restrictions  of  the 
deterministic  model  may  not  be  acceptable. 

5.  A  STOCHASTIC  ROUTE  FUNCTION  MODEL 

In  general,  a  stochastic  description  of  the  ensemble  of  routes  requires  the  joint  probability  density  on 
the  coordinate  positions  {9(x1,y1),...,9(x„,y,;)}  for  all  possible  samples  of  the  distances  traveled 
{x\  <  . . .  <  xn ;  7;  >  9} .  In  this  section,  we  present  a  stochastic  model  based  on  the  assumption  that  the  route 
set  probability  law  is  Markov.  We  first  specify  the  route  set  probability  law,  describe  how  the  probability 
law  is  used  to  generate  the  ship  tracks  in  terms  of  the  track  parameter  (90,y0,v),  and  conclude  by 
presenting  the  rate  function  that  specifies  the  Poisson  probability  law  on  the  ship-track  parameters. 

5.1  The  Route  Set  Probability  Law 

For  a  Markov  process,  the  joint  probability  density  is  completely  determined  by  the  “first-order” 
probability  density  and  the  probability  density  transition  function.  The  first-order  probability  density 
Pq(X]  ),y(x)  (9,  y; x)  describes  the  coordinates  of  all  routes  in  the  route  set  that  have  a  fixed  value  of  the 
distance  traveled  x.  In  general,  this  density  depends  on  x  (e.g.,  for  small  x  the  distribution  of  the  ship 
coordinates  is  concentrated  near  the  departure  port);  as  x  increases,  the  distribution  of  the  coordinates 
migrates  towards  the  destination  port  and  become  more  diffuse.  Note  that  at  the  departure  port  (x  =  9  )  the 
nominal  route  coordinate  9  is  zero,  so  that  the  first-order  density  has  the  form 


Pe(0),y(0)  (9,  y ;  x  =  0)  =  pJd  (yd  )5(0) . 


(26) 


In  the  following,  we  refer  to  y(j  and  Py  (y^  )  as  the  departure  cross-sectional  coordinate  and 
probability  density,  respectively. 

The  transition  probability  density,  denoted  by  Pq(x)  y(V)|0f  r')  y(x')  (®’ 9', y xc, xc') ,  is  the  probability 

density  on  the  coordinates  at  the  distance  x>  x  ,  (9(x),y(x)),  given  the  values  of  the  coordinates  at  the 
distance  x,  (O(x')  =  0',y(x')  =  y') .  The  first-order  density  for  the  coordinates  at  x  is  determined  from  the 
first-order  density  of  the  coordinates  at  x  by 


^e(x),y(x)  (0^*)  =  f  f  ^e(x),y(x)|e(x'),y(x')  (e,y;9',y',x,x>0(x'),y(x')  (e,,y,;x')rf0'rfy'  • 

0’  i 


(27) 


For  the  special  case  of  x  =  0 ,  the  transition  density  can  be  written  as 


770(x),y(x)|e(O),y(O)(0’7;^Y'  =  Yrf^^'=O)  =  ^0(x),y(x)lyrf(0’7;0>rf^)6(0')  ’ 


(28) 
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where  Pq^  y(x)\yd  (0S  Y5  0f,  >  *)  ’s  referred  to  here  as  the  departure  transition  probability  density. 

Substituting  from  the  last  equation  and  Eq.  (26)  into  Eq.  (27)  results  in 

/>0(*),y(x)(9.Y;*)=f  PQ(x),y(x)\ld{^r,Td^)Pyd{yd)dyd-  (29) 

y' 


5.2  The  Route  Set  Function 

As  with  the  deterministic  model,  the  route  and  the  track  functions  are  specified  in  terms  of  the  track 
parameter  (90,y0,v) .  The  set  of  real  ship  coordinates  SJ?,.  and  virtual  ship  coordinates  SRV  have  the  same 
form  as  for  the  deterministic  model,  where  the  route  envelopes  are  determined  by  the  support  of  the  first- 
order  probability  densities.  The  definition  of  the  route  function  is  as  follows. 

For  the  real  ship  coordinates,  the  route  function  is  determined  recursively  in  terms  of  the  transition 
probability  density.  The  process  is  as  follows.  At  time  t0 ,  the  initial  coordinates  are  specified  by  the  track 
parameter  (90,y0,v).  After  a  time  8 1 ,  the  ship  has  traveled  a  distance  8x  as  determined  by  the  ship- 
motion  function.  The  ship  coordinates  corresponding  to  this  increment  in  the  distance  traveled  are 
distributed  according  to  P0(Xo+5x)iy(Xo+5.Y)|e(Xo)iy(Xo)(9,Y;9'>Y,Ao  +  8x,x0) .  An  application  of  a  random 

number  generator  for  this  probability  density  yields  specific  values  of  the  ship  coordinates 
(9(x0  +  Sx),y(x0  +  Sx)).  The  process  is  repeated  to  yield  the  sequence  of  ship  coordinates  corresponding 
to  the  sequence  of  distances  traveled  determined  from  the  sequence  of  time  increments. 

For  the  virtual  coordinates,  we  take  y  =  yrf  and  9  =  x  in  analogy  with  the  deterministic  model.  At  the 
departure  port  (x  =  9  =  9 ),  the  ship  coordinates  are  (9,  y^) ,  so  that  the  recursive  process  described  above 
is  initiated  with  the  transition  density  Pq^  Y(x)|yrf  (9,y;y d',x)  •  As  such  the  cross-sectional  coordinate 

and  the  density  Py\yd(ld)  correspond  to  the  route  parameter  co  and  its  density  pm( co). 

5.3  The  Track  Parameters 


The  track  parameter  (90,y0,v)  is  related  to  the  track  parameter  (x0,yd,v)  =  (x0,eo,v)  by  the  first- 
order  probability  density  of  Eq.  (29).  It  follows  from  the  probability  law  of  Section  3  and  the  stochastic 
transformation  property  of  Appendix  A,  that  the  track  parameter  is  described  by  a  Poisson  process  (see 
Appendix  B).  For  the  real  ship  coordinates,  the  rate  function  is  given  by 


L(ld) 

fi0o,yo,v  (9o)Yo’vTo  )  —  J"  J"  Ax0  ,yd,v  (xo  T/d^v^o)  PQ(x),y(x)\  yd  (00’7o  ’xo’Yd  )dx0djd  ■>  (30) 

e/(  0)  0 


where 


^xo’Yd’v  iXo’Yd’  v’^o  )  fix  (jo  {xo’  Id  ’V’^o  )) 


>  Y  d  ATo  ) 


dxn 


(31) 


For  the  virtual  ship  coordinates, 

^o,,oAQo^o^t°,T)-pu  {y°)Pv{v)v  Px(^0-9ov  j,  (90,y0)esJlv(r). 


(32) 


18 


Richard  Heitmeyer 


Note  that  the  rate  function  of  Eq.  (31)  is  the  same  as  the  first  factor  of  Eq.  (18) 

with  co  replaced  by  and  Eq.  (32)  is  the  same  as  Eq.  (24)  with  |0 (y0 ;  0O  =  0)  replaced  by  pyd  (y0)  . 

The  application  of  these  equations  to  determine  the  shipping  densities  and  the  component  probability  laws 
is  the  same  as  that  described  in  the  preceding  section. 

6.  SUMMARY 

This  report  has  presented  two  models  that  describe  the  tracks  of  all  ships  en  route  in  a  region  during 
an  arbitrary  time  interval  [t0  ,t0+T].  Each  of  the  models  consists  of  the  composition  of  port-pair  models 
for  all  pairs  of  ports,  pseudo  as  well  as  real,  that  support  shipping  within  the  region.  Each  port-pair  model 
consists  of  a  track  function  that  describes  the  tracks  of  all  ships  present  during  [t0 ,ta+T]  and  a 
probability  law  on  the  number  of  ships  present  in  the  interval,  the  initial  positions  of  those  ships,  and  their 
nominal  speeds. 

The  probability  law  on  the  track  parameters  is  obtained  under  the  assumptions  that  the  times  at  which 
the  ships  depart  each  port  are  distributed  as  a  Poisson  process  with  a  time-varying  departure  rate  and  that 
the  ship  speeds  and  the  routes  that  the  ship  travel  are  statistically  independent.  Under  these  assumptions,  it 
is  shown  that  the  ship-track  parameters  are  distributed  as  a  space-time  Poisson  process  with  a  time- 
dependent  rate  function  determined  by  the  departure  rate  and  the  route  coordinate  density  functions.  This 
process  can  be  viewed  as  the  composition  of  two  Poisson  processes,  the  “en  route  process”  that  describes 
the  ships  that  are  en  route  at  t0  and  the  “entry  process”  that  describes  the  ships  that  depart  during  the 
interval  (t0,t0  +  T] .  The  en  route  process  is  defined  on  the  set  of  real  ship  coordinates;  the  entry  process 
is  defined  on  a  set  of  virtual  ship  coordinates.  The  rate  function  of  the  en  route  process  determines  a 
shipping  density  that  represents  the  mean  number  of  ships  per  unit  area  in  the  route  set  at  time  t0 .  A 
shipping  distribution  for  the  region  is  determined  by  integrating  the  shipping  density  over  the  resolution 
cells  in  a  longitude-latitude  grid.  For  both  models,  the  rate  function  specifying  the  Poisson  process 
depends  on  the  port-departure  rate  function  and  the  characteristic  speed  probability  density.  For  the 
deterministic  model,  the  rate  function  also  depends  on  the  route  coordinate  probability  density.  For  the 
stochastic  model,  the  rate  function  also  depends  on  the  transition  probability  density  of  the  Markov 
process. 

For  both  models,  the  track  function  is  determined  from  a  route  function  and  a  ship-motion  function. 
In  the  deterministic  model,  the  route  function  is  determined  from  the  route  coordinate  probability 
densities  in  such  a  way  that  each  route  is  uniquely  determined  from  any  point  on  that  route.  As  a 
consequence,  the  model  is  limited  to  route  sets  where  each  route  “progresses”  from  the  departure  port  to 
the  destination  port  without  doubling  back  on  itself  and  no  two  routes  can  cross  one  another.  The  route 
coordinate  probability  density  can  be  estimated  from  any  sample  of  ship  coordinates,  regardless  of 
whether  those  coordinates  are  organized  into  specific  routes. 

In  the  stochastic  model,  the  ship  tracks  are  obtained  as  realizations  of  the  Markov  process  on  the  ship 
route  coordinates.  The  process  is  recursive.  Starting  from  an  initial  position,  a  distance  traveled  is 
determined  from  the  motion  function  for  a  specific  time  increment;  the  ship  coordinates  are  then 
determined  from  the  initial  ship  coordinates  and  the  distance  traveled  using  the  transition  probability 
density  of  the  Markov  process.  The  process  is  repeated  to  determine  new  coordinates  by  replacing  the 
initial  ship  coordinates  with  the  current  ship  coordinates.  In  the  stochastic  model,  there  are  no  inherent 
constraints  on  the  route  crossings  as  in  the  deterministic  model.  Furthermore,  ship  location  data  used  to 
estimate  the  Markov  process  transition  probability  density  must  necessarily  be  organized  into  routes. 


Ship  location  data  that  is  not  organized  into  routes  cannot  be  used  to  estimate  the  Markov  process. 
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Appendix  A 

TRANSFORMATIONS  OF  POISSON  PROCESSES 


In  this  Appendix,  we  define  a  multidimensional,  nonhomogeneous,  Poisson  process  and  derive  two 
transformation  properties  that  are  used  in  the  development  of  the  probability  law  on  the  track  parameter. 
The  first  of  these  properties  applies  for  deterministic,  invertible  transformations;  the  second  applies  for 
stochastic  transformations. 

Al.  THE  POISSON  PROCESS 

Points  [zk  gZJ  are  distributed  as  a  Poisson  process  if,  for  any  two  disjoint  subsets  of  Z ,  Cj  and  C2 , 
the  number  of  points  in  Q  and  C2  are  statistically  independent  Poisson  random  variables  Nz  [Cj  ]  and 
Nz[C2\.  The  probability  mass  function  for  the  number  of  points  in  a  set  CcZ  is  given  by 

Pr{Az  [C]  =  n]  =  exp  {-AC  [C]}  — ,  (Al) 


where 


MZ[C]  =  \  [iz{z)dz  (A2) 

C 

is  the  mean  number  of  points  in  C  and  pz  (z  j  is  the  rate  function  for  the  process.  The  Poisson  process 
has  the  property  that  if  there  are  n  points  in  the  set  C  ,  then  the  positions  of  those  points  [zk;  k  =  1,...,«} 
are  independent,  identically  distributed  random  variables  with  common  probability  density, 

Pz\Nz[cynM  =  »M)Mz[CYX  •  (A3) 

Clearly,  to  specify  a  nonhomogeneous  Poisson  process,  it  suffices  to  specify  its  rate  function.  In  the 
special  case  where  the  rate  is  independent  of  z ,  the  Poisson  process  is  said  to  be  homogeneous.  In  this 

case,  Mz  [C]  =  pz  x  Volume [C]  and  pz  v  |y-]=,j  (zk)  =  Volume [C]  1 .  in  the  one-dimensional  case, 

where  Z  is  the  time  axis  and  C  is  a  time  interval,  C -\t0,t0  +T],  Eqs.  (Al),  (A2),  and  (A3)  take  the 
form  of  Eqs.  (10),  (1 1),  and  (12)  in  the  text. 

For  later  reference  we  note  that  the  characteristic  function  of  a  Poisson  process, 
cN,[C]  (“)  =  E\exp[iuNz  [C]}]  ,  is  given  by 

cN:[C]  (u)  -  cxp{Mz  [C](exp{i’w}  - 1)}  (A4) 
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and  that  a  stochastic  process  is  a  Poisson  process  if  its  characteristic  function  has  the  form  of  Eq.  (4)  in 
the  text. 

A2.  DETERMINISTIC  TRANSFORMATION  PROPERTY 

Let  the  points  in  the  space  Z  be  distributed  as  a  Poisson  process  with  density  function  pz  (z)  and  let 
g(z;g)  be  an  invertible  transformation  from  Z  to  a  space  Z'  for  each  value  of  the  parameter  q  .  Denote 

the  inverse  of  this  transformation  by  g  1  (z'lq)  .  Then  the  transformation  g  induces  a  Poisson  process  on 
Z'  with  a  rate  function  given  by 

Pz'(z')=ilz(?_1(z';^))*/(z';?)  >  (A5) 

where  j(z';g)  is  the  Jacobian  of  the  transformation. 

The  proof  of  the  transformation  property  is  as  follows.  Let  ScZ  be  any  set  and  let 
Aq  =  jz;  z  =  g-1  (z';g),  z'  g  c  Z  be  the  pre-image  of  B  under  the  transformation  g  .  Let  N'  be  the 
number  of  points  in  the  set  B  .  We  show  that  if  is  a  Poisson  random  variable  with  mean 

Mz'=\  Pz'(zVz'  >  (A6) 

B 

Where  pz-(z')  is  given  by  Eq.  (5).  It  suffices  to  show  that  the  characteristic  function  of  N' , 
cN'  ( u )  =  E [exp{mAz}J ,  can  be  written  in  the  form  of  Eq.  (4)  in  the  text.  To  this  end,  we  note  that  N' 
can  be  written  in  the  form 


W'=Z  Vz*) 

k 


(A7) 


where 


!Aq{zk) 


1  if  z  <=  A.. 

H 

0  otherwise 


is  the  indicator  function  of  the  set  Aq .  Substituting  from  Eq.  (A7)  into  the  definition  of  the  characteristic 
function  yields 
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Expanding  the  expected  value  operator  in  Eq.  (A8)  yields 


ctV'(u)  =  Z  Ezh...,zn\Nz 


exp  IAUk)\ 


Pr  {Nz=n}, 


(A9) 


where  Pr{Az  =  n  j  is  given  by  Eq.  (Al).  Since  the  points  {z/( ;  k  -  1 are  conditionally  independent 

and  identically  distributed  with  the  probability  density  of  Eq.  (A3),  the  expected  value  operator  inside  the 
summation  of  Eq.  (A9)  can  be  written  as 


Z]  >•  •  |  S 2  —n 


^  \  X 

k=n  ,  . 

exp U«Z  JA  {zk)\ 

Ezh...,zn\  Nz=n 

n  exP  1  iuI Aa  ( zk )  1 

V  k  \_ 

_k= 1  1  _ 

k=n 


=  n  Ezh..,zn\Nz=n  exP WlA  (zk) 
k= 1  L 


=  {E^n\Nz=n[Q^{MAq{zk) 

The  expected  value  in  the  last  equality  on  the  right-hand  side  of  Eq.  (A  10)  can  be  written  as 

E;kn\Nz=nWv 


(A10) 


iuIAq  {zk )}  =  Pr  {Aq  \NZ  =  n)  exp  [in]  +  (l  -  Pr  {Aq  \Nz  =  «}) 
=  Pr|^  | Nz  =  «j(exp{/w}  -l)  +  l. 


(All) 


Substituting  Eq.  (All)  into  the  last  equation  in  Eq.  (A10)  and  the  result  into  Eq.  (A9)  yields 

\n  exp 


ctV'(«)  =  Z  {^x{Aq\Nz  =«}(exp{m}-l)  +  l) 


:cxp. 


-Mjexpj^Prj^  | Nz  =  «j(exp{/w}  -l)  +  ljM 

or,  equivalently, 

cN>  (u)  =  exp  jdTPr  | Nz  -  wj  (exp {/w}  —  l) 

where  M  is  given  by  Eq.  (A2)  with  A  identified  with  Aq  .  Using  Eq.  (A3)  to  write 

■\Aq\Nz=n}=\  ^ dz 


(A  12) 


Pr 


(A  13) 
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and  substituting  into  the  right-hand  side  of  Eq.  (A12)  yields 


CN'  (w)  =  expj 


f  Mz) 


\Aq 


fife(exp{/«}  -l) 


(A  14) 


Equation  (A14)  indicates  that  N'  is  a  Poisson  process  with  mean 

Mz'=\  Mz)^z-  (A15) 

4 

Making  the  change  of  variables  in  the  integral  determined  by  the  transformation  g  yields 

Mz'  =  l  Vz(g~\z''’<lj)J{z'’<l)dz'  ,  (A  16) 

B 

as  was  to  be  shown. 

A3.  STOCHASTIC  TRANSFORMATION  PROPERTY 

Let  the  points  in  the  space  Z  be  distributed  as  a  Poisson  process  with  density  function  pz  (z) . 
Assume  that  given  any  sequence  of  points  [zk  e  Z;  k  =  1, ,  there  is  a  sequence  of  random  variables 
{zjeZ';i  =  l,...,n}  that  are  statistically  independent  and  identically  distributed  with  a  common 
probability  density  pzyz  [z'k',zk),  i.e., 

n 

^.i,...,4h,...,z„(zi'’---4;zi>---z«)  =  n  Pz'M^k)  ■  (A17) 

k= i 

Then  the  points  in  the  space  Z'  are  distributed  as  a  Poisson  process  with  a  density  p2-  (z')  given  by 

Pz'(z')=j  tv|z(z';z)Mz)Jz  •  (A18) 

z 

The  proof  of  the  transformation  property  is  as  follows.  Let  ScZ'  be  any  set  and  let  N'  be  the  number 
of  points  in  the  set  B  .  We  show  that  N'  is  a  Poisson  random  variable  with  mean 

Mz.  =  \  vAz')dz',  (A19) 

B 

where  py(z')  is  given  by  Eq.  (A18).  It  suffices  to  show  that  the  characteristic  function  of  A', 

cN'  («)  =  is  [exp{/wAf  }] ,  can  be  written  in  the  form  of  Eq.  (A4)  with  the  mean  given  by  Eqs.  (A  18)  and 

(A19).  To  proceed,  note  that  a  point  z'  is  in  the  set  B  if  there  is  a  point  zeZ  that  maps  into  the  point  z  . 
Thus,  the  random  variable  N'  can  be  written  in  the  form 
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tf'  =  X  IzxBizk> 4)  ,  (A20) 

k 

where 

fl  if  z  e  Z  and  z' e  i? 

Zx5(  k’  k)  j0 

is  the  indicator  function  of  the  set  ZxB .  But  Eq.  (20)  is  of  the  same  form  as  Eq.  (A7)  with  IA  ( Zj,  ] 
replaced  by  IZxB  { ) .  Thus,  proceeding  in  direct  analogy  to  Eqs.  (A8)  through  (A12),  we  obtain 

c^'(«)  =  exp|MPr|Zxi?|Az  =  «}(exp{/«} -l)j  ,  (A21) 

where  M  is  given  by  Eq.  (A2)  with  A  identified  with  Z .  Now, 

Pr{Zx5|Az  =«}  =  J  Pr{5|zeZ,Az  =n}p^Nx=n(z)dz  =  \  j  pz\N~=n{z)dz'  ^^-dz  ,  (A22) 

Z  “  Z  LB  "  J 

where  we  have  used  Eq.  (A3)  and  the  fact  that 

Pr{5|zeZ,Az  =«}  =  J  pz\Nz=n  (z)dz' .  (A23) 

B 

Substituting  from  Eq.  (A22)  into  Eq.  (A21),  we  obtain  Eq.  (A4)  with  pz  (z)  given  by  Eq.  (A19),  as 
was  to  be  shown. 


Appendix  B 

DERIVATION  OF  THE  TRACK-PARAMETER  PROBABILITY  LAWS 


In  this  Appendix,  we  derive  the  probability  law  on  the  track  parameter  (0o,yo,v)  for  both  the 
deterministic  and  the  stochastic  models.  The  probability  law  on  the  track  parameter  (A,0,(|)0,v)  is 
determined  from  the  probability  law  on  the  track  parameter  (0o,yo,v)  and  the  invertible  transformation 
relating  these  parameters.  To  obtain  these  results,  we  first  derive  the  probability  law  on  the  parameter 
( xQ ,  0),  v)  and  then  use  the  transformation  properties  of  Appendix  A.  Note  that  it  suffices  to  present  these 

derivations  for  the  single  port-pair  model  since  the  shipping  between  different  port-pairs  is  assumed  to  be 
statistically  independent.  We  conclude  this  Appendix  with  the  formal  definition  of  the  shipping 
distribution. 

Bl.  THE  TRACK-PARAMETER  PROBABILITY  LAW  FOR  (jc0,co,v) 

The  probability  law  on  the  track  parameter  (x0,co,v)  is  determined  from  the  probability  law  on  the 
track  parameter  (x,co,v)  and  the  transformation  defined  in  the  text.  To  derive  this  probability  law,  we 
first  show  that  the  points  |(x,co,v)|  are  Poisson  distributed  and  then  use  the  transformation  property  of 
Eq.  (A5)  in  Appendix  A  to  show  that  the  points  |(x0,co,v)|  are  Poisson  distributed  and  to  obtain  an 
expression  for  their  rate  function. 

To  this  end,  let  Z  be  the  set  of  pointsjz  =  (x,oc>,v);  x  e  (-00,00),  coe  [0,1],  v>  0| .  By  assumption,  the 
route  parameters  {<%}  and  the  ship  speeds  { v/(  j  are  statistically  independent  and  identically  distributed 
and  the  departure  times  j  x /. }  are  distributed  as  a  Poisson  process  with  rate  function  pT  (t ) .  Consequently, 
the  points  zeZ  are  distributed  as  a  Poisson  process  with  rate  function 

Mz)  =  Mx)Ad(®)/a(v)-  (B1) 

Next,  define  the  space  Z'  to  be  the  sets  of  points {z'  =  (x0,co,v)j  .  The  probability  law  on  Z'  is 
obtained  from  the  transformation  property  with  the  transformation  g(z,g)  of  Eq.  (A5)  determined  from 
Eqs.  (6)  and  (7)  in  the  text  and  the  identity  mappings  on  CO  and  v  .  Specifically,  g(z,q )  is  determined  by 

x0  =  Mx(t0-x;m,v)  ,  (B2) 

where  the  initial  time  t0  in  Eq.  (B2)  is  the  parameter  q  of  the  transformation  g(-,q).  The  inverse 
transformation  g_1(z',g)  is  determined  by 

*=*<,  (B3) 
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The  Jacobian  of  the  transformation  is 

j(xo,(0,v;to)- 


5?o  (  Xq  ’  l0  ) 


5x„ 


Applying  the  transformation  property,  we  conclude  that  the  points  |z' =  (x0,eo,v)  e  Z'j  are  Poisson 
distributed  with  rate  function  given  by 

Px0,co,v (^■05®,v>fo,7’)  —  Ax0,oo,v (co)/?v  (v)  ,  (B4) 

where 

I  (•*’0’®'  A  ) 


A.YQ  ,CO,V  (*0  J  A  ^)- 


Pt(^0  *30  (^0>  ®),V,f0  )) 


OX, 


for  0<xo  <L(co) 


for  -vT  <x0<  0 


(B5) 


B2.  THE  TRACK  PARAMETER  PROBABILITY  LAW:  DETERMINISTIC  MODEL 

The  probability  law  on  the  track  parameter  (0o,yo,v)  is  determined  from  the  probability  law  on  the 
track  parameter  ( x0 ,  co,  v)  and  the  transformations  defined  in  the  text.  To  derive  this  probability  law, 
define  the  spaces  Z'  and  Z"  to  be  the  sets  of  points  |z'  =  (x0,  oo,v)|  and  |z"  =  (0o,yo,v)} ,  respectively. 
The  transfonnation  from  Z'  to  Z"  is  obtained  from  Eqs.  (4)  and  (15)  in  the  text  with  x  replaced  by  x0  ; 

0O  —  0(xo,co) 

Y0  —  (ro,0(xo,co)) 

The  inverse  transformation  is  given  by  Eqs.  (16)  and  (17)  in  the  text: 

®  =  py|e(Y0;90) 
x0  —  X 0  (0o,yo )  -  ^(00,®) 


(B6) 


(B7) 


<i>=P. 


|e(yo;0o) 


The  Jacobian  of  the  transformation  is 


j(xo,(0,v,to)  = 


dX0{Q0,y0)  dX0  (0o,yo ) 


50, 


sy c 


DR 


|0(yo;0o)  5Py|0(yo;0o) 


50, 


sy c 


To  evaluate  the  Jacobian,  we  expand  the  determinant  and  then  use  the  chain  rule  to  evaluate  the 
partial  derivatives  of  Xo(0o,yo),  where  co  is  given  by  Eq.  (13)  in  the  text.  Finally,  we  use  Eq.  (B7)  to 
obtain 


J  ( x0 ,  oi,  v,  t0 )  p(Yo»90) 


dX0{Q0,y0) 


50, 


(B8) 
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It  follows  from  the  transformation  property  of  Appendix  A,  Eqs.  (B4)  and  (B5),  and  the  fact  that  co  is 
uniformly  distributed,  that  the  points  \z"  =  (0o,yo,v)  e  Z" j  are  Poisson  distributed  with  the  rate  function 

given  by  Eq.  (18)  in  the  text.  Furthermore,  since  the  points  z"  are  distributed  as  a  Poisson  process,  the 
initial  ship  positions  (0o,yo)  are  distributed  as  a  Poisson  process  with  rate  function  given  by  Eq.  (19). 

Equation  (20)  in  the  text  follows  from  the  general  equation  for  the  mean  number  of  points  in  a  set  (Eq. 
(Al)  in  Appendix  A)  with  the  rate  function  given  by  Eq.  (19)  in  the  text.  Eq.  (21)  in  the  text  follows  from 
Eq.  (A3)  in  Appendix  A. 

To  obtain  the  probability  law  for  the  virtual  ship  coordinates,  we  note  that  (a)  these  ships  travel  with 
constant  speed  motion;  (b) X0 (0o,yo ) = 0O ,  and  p(y0',Q0 )  =  p(j0',Q0  =  0) ;  and  (c)  the  condition 

x  e  (t0 , t0  +  r)  is  equivalent  to  the  condition  x0  -  90  e  \-vT, 0) .  The  results  follow  using  the  same  logic 
as  that  used  to  obtain  the  probability  law  on  the  real  ship  coordinates. 

B3.  THE  TRACK-PARAMETER  PROBABILITY  LAW:  STOCHASTIC  MODEL 

The  probability  law  for  the  stochastic  model  is  determined  from  the  probability  law  on  (x0,oo,v)  and 
the  stochastic  transformation  property  of  Eq.  (A18).  This  follows  from  the  fact  that  the  track  parameter 
for  the  stochastic  model  is  just  (x0,co,v)  with  co  replaced  by  yrf  .  Thus,  (x0,y^,v)  is  Poisson- 

distributed  with  rate  function  given  by  Eqs.  (B4)  and  (B5)  with  co  replaced  by  y  j  and  pm  ( co)  replaced 

by  Py  (y0) .  To  apply  the  Stochastic  Transformation  property,  note  that  the  probability  density  in  Eq. 

(A18)  is  the  transition  probability  density  of  Eq.  (29)  in  the  text.  Substitution  into  Eq.  (A18)  completes 
the  derivation. 


B4.  THE  SHIPPING  DISTRIBUTION  FUNCTION 


The  shipping  distribution  represents  the  mean  number  of  ships  in  each  resolution  cell  determined  by  a 
rectangular  grid  defined  in  the  (A-,  4>)  coordinate  system.  The  formal  definition  of  this  function  is  as 

follows.  Let  ,  <J)y  be  the  rectangular  grid  of  points  determined  by  Xt  =  Xc  +  AX(i  - 1);  i  - 1,...,/ 

and  <\>j  =<|)c  +  Ac)) ( y  —  1 ) ;  i  -  1,...,  Js ,  where  AX  and  Acj)  are  the  longitude  and  latitude  spacings, 
respectively,  and  Xc  and  <t>c  are  the  coordinates  of  the  southwest  comer  of  the  region.  Further,  let 


C 


,  AX  .  AX 

Xi - ,X;+ - 

1  2  1  2 


A4>  Ac)) 

■tv-f-W 


(B9) 


be  the  resolution  cell  corresponding  to  the  point  j  and  let  C'[Xi,§  j  j  be  the  pre-image  of  that 

resolution  cell  in  the  auxiliary  coordinates  obtained  under  the  transformation  between  the  longitude- 
latitude  and  the  auxiliary  coordinates.  Finally,  let  Dn  m  (Xi,^J  )  be  the  shipping  distribution  for  the  single 

port  pair Pn  m  [Xpfyj ) .  Then  it  follows  from  Eq.  (B1 1)  that  Dn  m  (\-,<t>y )  is  given  by 


D 


n,m 


j  )  (j,  \P0,y ;n,m  {^o’To’^o  ^^0^1  c 


(BIO) 


The  shipping  distribution  for  the  total  region  follows  from  the  fact  that  an  aggregate  of  independent 
Poisson  processes  is  itself  a  Poisson  process  with  a  rate  function  and  shipping  density  given  by  the  sum  of 
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the  rate  functions  and  the  shipping  densities  of  the  component  processes.  Thus,  the  shipping  density  for 
the  region  is  given  by 

NP 

To,y  (®o’Yo’^o  )  —  ^Q,y,n,m  Yo 

n,m= 1 

and,  hence,  the  shipping  distribution  for  the  region  is  given  by 

N p 

Z  Dn,m(X^j\  (B12) 

n,m=\ 


Appendix  C 

A  SIMPLIFIED  DETERMINISTIC  MODEL 


As  seen  in  the  text,  the  probability  law  on  the  track  parameter  for  the  deterministic  model  is 
considerably  simplified  by  assuming  constant  speed  motion  and  by  approximating  dx/dQ0  by  unity.  In 
this  Appendix,  we  describe  further  simplifications  in  the  model  obtained  by  making  four  additional 
assumptions.  The  first  states  that  the  departure  rate  is  independent  of  time  (i.e.,  pT  (?)  =  pT ).  This 

assumption  is  often  necessitated  by  practical  considerations,  since  the  available  departure-time  data  may 
be  sufficient  to  estimate  the  time-averaged  departure  rate  but  not  its  time  dependence.  The  second 
assumption  states  that  the  route  coordinate  probability  density  /;.,  (y;9)  has  the  same  form  for  all  cuts  in 

the  route  set.  This  assumption  results  in  a  simplified  track  function.  The  third  assumption  provides  a 
simple  approximation  for  the  shipping  distribution  and  the  track-parameter  probability  density  for  a 
sufficiently  fine  grid  resolution  grid.  Finally,  the  fourth  assumption  states  that  the  route  set  itself  can  be 
approximated  as  a  sequence  of  simple  route  sets,  where  the  nominal  route  for  each  is  a  segment  of  a  great 
circle  arc.  This  assumption  results  in  simple  expressions  for  the  transformation  between  the  auxiliary 
coordinates  and  the  longitude- latitude  coordinates. 

Cl.  CONSTANT  DEPARTURE-RATE  PROBABILITY  LAW 

The  simplified  track-parameter  probability  law  for  the  real  ship  coordinates  is  obtained  by 
substituting  pT  for  the  time-dependent  departure  rate  in  the  expressions  for  the  track-parameter  rate 
function  and  the  shipping  density  (Eqs.  (22)  and  (23)  in  the  text)  and  then  substituting  the  results  into  the 
expression  for  the  mean  (Eq.  (20  in  the  text))  and  the  track-parameter  probability  density  (Eq.  (21}  in  the 
text)).  The  shipping  density  becomes 

M-0,y  (®o’Yo  )  —  .Py|0  (To’^o  )PtP  ’  (Cl) 

where 


P  =  j  v  lPv(v)dv  (C2) 

0 

is  the  mean  reciprocal  ship  speed.  Furthermore,  the  mean  number  of  ships  with  initial  positions  in  a  set 
C,=[0i,02]x[yi,Y2]  becomes 

_02 

M,.(c,.)  =  pTp{ 

01 

Note  that  the  mean  of  all  ships  that  are  en  route  at  time  t0  is  simply  Mr  (ill,. )  =  pTPT  .  Finally,  the  track- 
parameter  probability  density  becomes 

Pe  (0q ’  Yo  C,- )  —  pv  (v).Py|Q  (Yo ’0<3  5  (C4) 


V*2;0)-Vyi;0)  de 


(C3) 
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where 

A-(v)  =  (Pv)  '  Pv{v)  .  (C5) 

According  to  this  probability  density,  the  Q0  coordinate  is  uniformly  distributed  on  the  interval  [ 0 ,  A  ] 
and,  given  90,  the  y0  coordinate  has  the  probability  density  q(yo;0o)-  The  ship  speed  has  the 
probability  density  pv  ( v)  and  is  statistically  independent  of  the  ship  positions. 


For  the  virtual  ship  coordinates,  there  is  no  appreciable  simplification  in  the  shipping  density  other 
than  that  obtained  by  replacing  pT|t0  -0ov-1  j  by  pT  in  Eqs.  (24)  and  (25)  in  the  text.  Flowever,  for  sets 

of  the  form  Cv  -  [-oo,0]x[y1,y2],  the  mean  number  of  ships  simplifies  to1 


Me(T-Cv)  =  \ixT 


Py|e(Y2;90  -°)_-Py|e(Yf  0O  -°) 


(C6) 


Note  that  the  mean  number  of  all  ships  that  depart  in  the  interval  is  simply  Me  (T\^\v)  =  pxT ,  as 
expected  since,  by  assumption,  the  departure  times  are  Poisson  distributed  with  rate  pT.  Finally,  the 
track-parameter  probability  density  for  the  ships  that  depart  in  the  interval  (t0,t0+T~\  becomes 


Pd  ( j  Yo  ATo ) 


(Tv)  1  Pv(v)py^(y0;B0  =0)/L  for 0O  e [-Tv, 0] 
0  otherwise 


(C7) 


According  to  this  probability  density,  the  y0  coordinate  has  the  probability  density  (yo;0o  =  0) 

and  is  independent  of  both  the  90  coordinate  and  the  ship  speed.  The  90  coordinate  and  the  ship  speed 
are  statistically  independent  with  the  90  coordinate  uniformly  distributed  on  the  interval  [-vT,  0) .  The 
ship  speed  has  the  density  pv  (v) ,  rather  than  pv(v). 


Equations  (C3)  through  (C7)  specify  the  track-parameter  probability  law  in  the  simplified  model.  For 
those  ships  en  route  at  time  t0 ,  the  probability  that  there  are  n  ships  is  given  by  Eq.  (10)  in  the  text  using 
the  mean  of  Eq.  (C3)  and  the  track  parameters  for  those  ships  are  independent  and  identically  distributed 
with  the  common  probability  density  of  Eq.  (C4).  For  those  ships  that  depart  during  the  interval 

(t0,t0  +r] ,  the  probability  that  there  are  n  ships  is  also  given  by  Eq.  (10)  using  the  mean  of  Eq.  (C6), 

and  the  track  parameters  for  those  ships  are  independent  and  identically  distributed  with  the  common 
probability  density  of  Eq.  (C7). 

Specific  ship  track  realizations  can  be  obtained  using  the  track  function  with  specific  realization  of 
the  track  parameter  obtained  from  the  probability  law  of  Eqs.  (C3)  through  (C6).  In  particular,  for  those 
ships  en  route  at  time  t0 ,  a  Poisson  random  number  generator  can  be  used  to  determine  an  integer  n  that 
represents  the  number  of  ships  en  route  at  time  t0 .  Then,  for  each  of  these  n  ships,  random  number 
generators  for  the  three  probability  densities  of  (C4)  can  be  used  to  determine  the  corresponding  track 
parameters.  A  similar  procedure  can  be  used  for  the  ships  that  depart  in  the  interval  (t0,t0  +  T  1 .  An 

alternate  procedure  for  generating  track-parameter  realizations  is  obtained  under  the  fine  resolution 
assumption  described  below. 


1  The  integral  in  the  expression  for  the  mean  is  evaluated  by  interchanging  the  order  of  integration. 
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Finally,  the  shipping  density  of  Eq.  (Cl)  determines  the  shipping  distribution  for  the  longitude- 
latitude  grid  of  interest.  For  fine-resolution  grids,  where  the  cell  size  is  small  with  respect  to  the  variation 
in  the  shipping  density,  the  shipping  distribution  is  approximately  given  by  Eq.  (Cl 8). 


C2.  THE  INVARIANT  ROUTE  COORDINATE  PROBABILITY  DENSITY 


The  invariant  route  coordinate  density  assumption  states  that  ^7jg(yo;0o)  has  the  same  form  for  all 
values  of  90  .  This  is  equivalent  to  the  assumption  that  there  is  a  “normalized”  cross-sectional  variable  y 
and  a  probability  density  p0  (y) ,  such  that,  for  each  90,  the  probability  density  ( y0 ; 0o )  is 

determined  from  pa  (y)  by  the  transformation 


y  =  Ae(90  )y  +  ea  (0O );  for  ye  [-1,1] 


where 


Ae  =  (^(0o)-e/(9o))/2;  =(^(0o)+e/(9o))/2  • 

Under  this  assumption,  the  components  of  the  track  function  are  given  by2 
GQ(t-t0-,Q0-,y0-,v)=Q(t-t0)=v(t-t0)+%0 
Gy  (t  - 10\ 90 ;y0 ; v)  -  \ (y0  - e*  (90  ))c(90 )_1  A e(&(t -t0))  +  ea  (0(t - 10 )) 


(C8) 

(C9) 


(CIO) 


Note  that  this  track  function  is  independent  of  the  form  of  the  route  coordinate  density  function. 


To  establish  Eq.  (C10),  we  note  that  it  follows  from  the  transformation  of  Eq.  (C8)  that  the  route  set 
density  and  the  route  set  distribution  function  are  given  by 


Py\<d{r,%)  =  P0 

Py\ q{t^o)  =  Po 


Y~gfl(9) 
Aea  (9) 

Y~gfl(9) 

A  ea(Q) 


A^  (9)_1 


(Cll) 


where  p0  (y)  and  P0  (y)  are  the  probability  density  and  the  distribution  function  for  the  normalized 

cross-sectional  variable,  respectively.  Furthermore,  it  follows  from  the  second  equation  of  Eq.  (Cll)  and 
the  definition  of  the  cross-sectional  route  function  in  the  text  that,  for  each  route  eo ,  we  must  have 

[(y-efl(0))/Ae(0)]  =  Po-1(co)  ,  (C12) 


where  P0  1  is  the  inverse  of  the  distribution  function  P0  .  Now,  Eq.  (Cl 2)  must  hold  for  all  points  on  the 
route,  including  the  initial  point  (90,y0 ) .  Thus,  we  also  have 

[(T-gfl(0o))/Ag(0o)]  =  V(®)  -  (Cl 3) 

Comparing  these  last  two  equations  yields 


2  To  obtain  this  result,  it  is  necessary  to  identify  x  with#.  Although  this  follows  from  the  assumption  that 
d'Xj  dd„  =  1 ,  it  is  clearly  an  approximation  for  any  nondegenerate  route  set. 
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Y  =  (y0  - U,  (0O  ))(Ae(0)/Ae(0o ))  +  ea  (0) .  (C14) 

To  obtain  the  track  function,  we  identify  x  with  0  and  use  the  assumption  of  constant  speed  motion 
to  write 

0(t-to)  =  v(t-to)  +  0o  .  (Cl 5) 

Equation  (CIO)  follows  immediately  from  Eqs.  (Cl 4)  and  (Cl 5). 

We  conclude  this  subsection  by  noting  that  the  track  function  of  Eq.  (CIO)  can  be  further  simplified 
by  taking  the  lower  route  envelope  to  be  the  negative  of  the  upper  route  envelope,  i.e., 
e/( 0)  =  -ea(0)  =  -e(0).  For  this  case,  ea(0)  =  0  and  Ae(0)  =  e(0)  so  that  y  component  of  the  track 
function  of  Eq.  (CIO)  becomes 

Gy(t-to;0o;yo;v)  =  Y  =  [yo/e(0o)]Ae(0(t-to)).  (Cl  6) 


A  particularly  useful  choice  for  the  normalized  probability  density  p0  (y)  is  the  beta  density 
modified  to  span  the  interval  [-1, 1] ,  i.e., 


Paid) 


T(a  +  P) 

n  +  y) 

a-1 

( 1  -  yP 

2T(a)T(p) 

l  2  J 

l  2  J 

(Cl  7) 


The  route  coordinate  probability  density  used  for  the  San  Diego  ship  track  example  in  Ref.  5  is  the 
Beta  density  for  a  =  [3  =  2.  The  route  envelope  u(0)  was  taken  to  be  a  segment  of  a  sinusoid  with  an 
amplitude  that  depends  on  the  length  of  the  route  [6]. 


C3.  THE  FINE  RESOLUTION  ASSUMPTION 


We  assume  that  the  resolution  cells  c(Xj,$j )  in  the  grid  T  are  sufficiently  small  that  the  shipping 

density  is  approximately  constant  over  each  resolution  cell  and  equal  to  its  value  in  the  center  of  the 
resolution  cell.  Under  this  assumption,  the  port-pair  shipping  distribution  is  approximately  given  by 


D 


(h’$j  )  =  E0,y  (®(X,f ,  <t>y  )  ,T  (*,,■  ,§j  ))  A  (<t>y  )  , 


(CIS) 


where  is  the  area  of  the  resolution  cell  C[k-n§j  j  and  ©(z^,^  j  and  r  (z.^p^  j  are  the 

coordinate  functions  in  the  mapping  from  the  latitude-longitude  coordinates  to  the  auxiliary  coordinates. 
Substituting  for  the  shipping  density  from  Eq.  (Cl)  into  Eq.  (Cl 8)  yields 

D(h> <tV )  =  PxP^Y|e  (r (?q- , <t>y ) ; @ (X,,- , <|)y )) (<|)y )  •  (Cl 9) 

For  small  resolution  cells,  A(<\ )  A  is  approximately  given  by 

A($j)  =  Re2  cos  ((f)y )  AkA(f)  (C20) 


where  Re  is  the  radius  of  the  Earth. 


To  obtain  an  alternate  means  of  generating  the  ship-track  parameters,  we  derive  the  probability 
density  on  the  track  parameter  given  that  the  ship  is  in  the  resolution  cell.  This  probability  density  is  given 
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by  the  rate  function  divided  by  the  mean  number  of  ships  in  the  cell.  Under  the  fine  resolution 
assumption,  the  mean  Mr  ( tQ;C )  is  given  by  Eq.  (Cl 8)  and  the  rate  function  is  constant  in  the  resolution 

cell  with  the  value  determined  by  replacing  90  and  y0  by  )  and  r(\-,cj)y- ) ,  respectively.  Thus, 

we  have 


Pe(h’$j’V’C) 


Vk&v  (0  )  ,■ r  Aj )  .v) 


(C21) 


For  the  simplified  model  shipping  density  of  Eq.  (Cl),  this  probability  density  becomes 

£v(v) 


Pe(h’§j’v’C)- 


A 


(O  ’ 


(C22) 


where  pv  (v)  is  given  by  Eq.  (C5). 

According  to  Eq.  (C22),  given  that  a  ship  is  in  the  resolution  cell,  its  position  in  that  cell  is  uniformly 
distributed  and  its  speed  is  independent  of  its  position  and  has  the  probability  density  pv  (v) .  Thus,  under 

the  fine  resolution  approximation,  the  ship  track  realizations  can  also  be  generated  by  first  using  a  Poisson 
random  number  generator  with  the  mean  given  by  Eq.  (C6)  to  determine  the  number  of  ships  in  the 
resolution  cell  and  then  using  uniform  random  number  generators  for  the  probability  densities  of  Eq. 
(C22)  to  generate  the  track  parameters. 

C4.  COMPOUND  ROUTE  APPROXIMATION 

The  transformation  between  the  auxiliary  coordinates  and  the  longitude-latitude  coordinates  has  a 
particularly  simple  form  in  the  special  case  where  the  nominal  route  is  the  great  circle  arc  connecting  the 
departure  port  with  the  destination  port.  Elowever,  for  many  port-pairs,  it  is  not  possible  to  choose  the 
nominal  route  to  be  a  great  circle  arc  since  the  routes  must  bend  to  avoid  land  masses.  For  these  port- 
pairs,  the  nominal  route  can  be  approximated  as  a  sequence  of  great  circle  arcs  connected  at  their 
endpoints  (see  Fig.  Cl).  For  each  of  these  great  circle  arcs,  the  routes  can  be  described  in  the  auxiliary 
coordinates  using  a  route  function  defined  for  the  great  circle  arc  and  the  simple  coordinate 
transformation  can  be  used  to  map  those  routes  to  the  longitude-latitude  coordinates.  These  “simple” 
route  sets  can  then  be  joined  together  to  form  the  complete  route  set  between  the  departure  port  with  the 
destination  port.  A  method  for  connecting  simple  route  sets  to  form  a  “compound  route  set”  is  described 
in  Springer  et  al.3  along  with  the  coordinate  transformations  for  the  great  circle  nominal  routes. 


P.  Scrimger,  R.  Heitmeyer,  and  P.  Boulon,  “A  Computer  Model  of  Merchant  Shipping  in  the  Mediterranean  Sea,” 
SACLANT  Undersea  Research  Center  Report,  SR- 164,  La  Spezia,  Italy,  SACLANT  Undersea  Research  Center, 
1990. 


Richard  Heitmeyer 


